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The surface critical behavior of semi-infinite (a) binary al- 
loys with a continuous order-disorder transition and (b) Ising 
antiferromagnets in the presence of a magnetic field is con- 
sidered. In contrast to ferromagnets, the surface universality 
class of these systems depends on the orientation of the sur- 
face with respect to the crystal axes. There is ordinary and 
extraordinary surface critical behavior for orientations that 
preserve and break the two-sublattice symmetry, respectively. 
This is confirmed by transfer-matrix calculations for the two- 
dimensional antiferromagnet and other evidence. 
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A cornerstone of the modern theory of critical phenom- 
ena is the idea of distinct universality classes jjj. The 
bulk universality class to which a particular system be- 
longs is generally determined by a few basic properties, 
such as the spatial dimension and the symmetries of the 
order-parameter. In the renormalization-group picture 
each universality class corresponds to the basin of at- 
traction of a fixed point in a many-dimensional space of 
Hamiltonians. 

Over the past two decades convincing evidence has 
emerged for similar universality classes in the surface crit- 
ical behavior of semi-infinite systems close to the bulk 
critical point |@,0,H|. The surface universality class of a 
particular system is determined by (i) the bulk univer- 
sality class and (ii) additional relevant surface properties. 
In a ferromagnet in zero magnetic field, for example, the 
strength of the spin couplings near the surface is one 
of these additional properties. For subcritically, crit- 
ically, or supercritically enhanced surface interactions, 
the surface critical behavior is "ordinary" , "special" , or 
"extraordinary", respectively. Another relevant surface 
property is a surface ordering field. The surface criti- 
cal behavior of a ferromagnet with zero bulk magnetic 
field and nonzero surface field belongs to the "normal" 
universality class || , independent of the strength of the 
surface couplings. In both the extraordinary and normal 
cases there is long-range order near the surface above the 
bulk critical temperature, and both transitions have the 
same universal properties f|. Following common usage 
we refer to the joint universality class as extraordinary. 
Note that the Ising model with d < 2 does not exhibit 
a "true" extraordinary transition (except for infinite sur- 
face couplings), because of the low boundary dimension, 



but there is a normal transition. 

Field-theoretic studies of Ising surface critical behavior 
p| begin with the one-component </> 4 model with Hamil- 
tonian 

defined on the half space ¥1+ = {(x\\,y) e M d | y > 0} 
with boundary plane S at y = 0. In this model there 
is ordinary, special, and extraordinary, surface critical 
behavior at bulk criticality for hi = and cq > c sp , 
Co = c sp , and cq < c sp respectively. For h\ ^ and 
arbitrary cq the surface critical behavior is normal . 

In this Letter the surface critical behavior of (a) binary 
alloys with a continuous order-disorder transition and of 
(b) Ising antiferromagnets in the presence of a magnetic 
field H is considered. Our main result is that these sys- 
tems differ from ferromagnets in an important respect. 
The surface universality class depends on the orientation 
of the surface plane with respect to the crystal axes. 

That the orientation affects the surface behavior of sys- 
tems (a), (b) has been recognized by Schmid Q, who 
carried out Monte Carlo simulations and mean-field cal- 
culations for a lattice model of the A2-B2 order-disorder 
transition j9) in FeAl, equivalent to an Ising antiferro- 
magnet with a free surface on a bcc lattice. According to 
Ref. (^] , for nonideal bulk stoichiometry (nonzero mag- 
netic field in the equivalent Ising antiferromagnet) the 
order-parameter symmetry is broken by the (100) sur- 
face orientation, and long-range order persists near the 
surface above the bulk critical temperature T c . For the 
(100) surface and ideal stoichiometry or the (110) surface 
and arbitrary stoichiometry, the symmetry is not broken, 
and both the bulk and surface order vanish for T > T c . 
We agree with these results but question the conclusion 
that the surface critical behavior at T c is ordinary in all 
of the above cases. We find that the surface universality 
class depends on the orientation of the surface, with ex- 
traordinary surface critical behavior for orientations that 
break the symmetry of the A and B sublattices (consis- 
tent with the renormalization-group relevance of surface 
ordering fields) and ordinary surface critical behavior for 
symmetry-preserving orientations. Transfer-matrix cal- 
culations in d — 2 described below confirm these predic- 
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tions. The predictions are also consistent with an effec- 
tive field-theoretic Hamiltonian derived from mean-field 
theory, as outlined below. 

As in [H, we consider a simple lattice-gas model in 
which each site of a bcc lattice is occupied by a particle of 
type 1 or 2. Nearest-neighbor pairs contribute Vn, V22, 
or V12, depending on the species involved, to the total 
energy. In the grand canonical ensemble with chemical 
potentials /l«i,/X2 the lattice gas is equivalent to an Ising 
model with Boltzmann factor exp(— 7ii a t), where 
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The bulk and excess surface magnetic fields are given by 
H = - M2) - K^n " V22)](k B T)- 1 and H x = 

j((-(i)(Vu-V 22 )}(k B T)~\ where ( and Ci are the co- 
ordination numbers of interior and surface sites, respec- 
tively. The nearest-neighbor coupling K = \{Vu + V22 — 
2Vi2)(A:_Br) _1 is assumed to be positive, corresponding 
to an antiferromagnet, and is the same for surface and 
interior spins. For \H\K~ 1 < C, \H + H^K' 1 < &, 
the system is antiferromagnetically ordered at T = 0. 
The two-dimensional analog of the model is shown in 
Fig. 1, with the A and B sublattices indicated by filled 
and empty points. 

+ 
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FIG. 1. Symmetry breaking (10) surface (vertical broken 
line) and symmetry preserving (11) surface (diagonal broken 
line) . 

The standard definition of the bulk order parameter of 
the antiferromagnet in a magnetic field is the difference 
of the sublattice magnetizations. The conjugate ordering 
field is a staggered magnetic field. The bulk field H in (||) 
is an example of a nonordering field. Increasing H low- 
ers the critical temperature, at least for weak fields, but 
does not alter the bulk universality class. As can be seen 
by redefining the spins on one sublattice with a minus 
sign, the Ising antiferromagnet is exactly equivalent to 
an Ising ferromagnet with a nonordering staggered field, 
irrelevant in the bulk. Thus the same universality classes 
as in the ferromagnetic Ising model are expected. 

The relevance of the surface orientation can be under- 
stood from Fig. 1 . Suppose that there are N — > 00 layers 
of spins, with free surfaces at the first and Ath layers and 
periodic boundary conditions in the other direction. For 
the (11) surface orientation adjacent spins in any layer 



belong to different sublattices. Both the bulk and sur- 
face fields H and H\ in Eq. (0) are nonordering. The 
Hamiltonian is invariant under a translation of all the 
spins parallel to the surface by one lattice constant, i.e., 
under interchange of the A and B lattices. Thus the (11) 
orientation respects the symmetry of the two sublattices. 
The order-parameter profile (f> n , defined as the difference 
of the sublattice magnetizations in the nth layer, van- 
ishes for T > T c due to symmetry of the sublattices. For 
T < T c there is ordering at the surface driven by the 
bulk order. Thus ordinary surface critical behavior is 
expected. 

A (10) surface breaks the A-B symmetry of the antifer- 
romagnet, since all the surface spins belong to a particu- 
lar sublattice. The order-parameter profile is symmetric 
and antisymmetric about the midpoint of a strip with N 
layers for N odd and even, respectively (see Fig. 3 and 
Fig. 4 below). This can be understood in terms of the 
equivalent ferromagnet. The bulk staggered field is con- 
stant within a layer but alternates in sign from layer to 
layer. The fields of magnitude \H + H\\ on the two sur- 
faces, which are parallel for N odd and antiparallel for TV 
even, are ordering fields for this orientation. Under the 
renormalization group the bulk staggered field is driven 
to zero, but surface fields with H — h or H — orientations 
survive. For both the H — h and H — orientation extraordi- 
nary critical behavior is expected in the limit N — > 00, 
due to the field-induced order at the surface. 

Applying the same reasoning to the three-dimensional 
analog of the above system (bcc lattice, N — > 00 layers, 
two free surfaces, and periodic boundaries in all other 
directions), we predict ordinary and extraordinary sur- 
face critical behavior for (110) and (100) surfaces, re- 
spectively. 

The above conclusions also apply if either of the two 
fields H, Hi vanishes. When H = Hi = (ideal sto- 
ichiometry) the magnetization vanishes on both sublat- 
tices for T > T c . Ordinary critical behavior is expected 
for arbitrary surface orientations. 

We have checked our predictions in d = 2 with numeri- 
cal transfer-matrix methods. First the correlation length 
£ N (K c (H),H,Hi) parallel to (11) and (10) surfaces of 
the antiferromagnet (Q) defined on a strip of square lat- 
tice with N layers of spins was calculated. Here K C (H) 
is the bulk critical coupling constant, determined as a 
function of H with great precision in Ref. [[l0| . Then the 
universality class was deduced by comparing irN^ 1 ^, 
with £tv expressed in units of the layer separation, with 
the known amplitudes 



A ab = iT lim N- 1 ^ 
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for ferromagnetic Ising strips with boundary conditions 
a, b. The A ab , which are clearly scale invariant, are uni- 
versal quantities that depend on the surface universality 
class. For ferromagnetic Ising strips with ordinary (zero- 
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field) boundary conditions and with parallel and antipar- 
allel surface fields 



A° 



= 2 



A ~ 2 



A + - = 1 
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The first two amplitudes in Eq. (Q) follow from the ex- 
act [|ll] results r;° rd = 1, ?7| x = 4 for the semi-infinite 
two-dimensional Ising model and Cardy's Q] formula 
A aa = 2(?7") _1 for homogeneous critical systems with 
the same boundary condition a on both edges. This re- 
lation is a consequence of conformal invariance and the 
conformal mapping of the half plane with boundary con- 
dition a onto the strip. Here the surface critical exponent 
?7|| is defined by the r~ v ^ decay of order-parameter cor- 
relations parallel to the boundary of the half space. The 
entry A + ~ = 1 in Eq. (jj) is derived in Ref. |Q. It also 
follows from the exact spin-spin correlation function |l3| ] 
for -| — boundary conditions. 

We have numerically analyzed the transfer matrices of 
rectangular slices of the strip with dimensions (N— 1) x r, 
where r = 1 for (11) edges and r = 2 for (10) edges. The 
correlation length £jv> expressed, like the width N — 1 
of the strip and r, in units of the layer separation, was 
determined from the largest and next largest eigenvalues 



of the transfer matrix \$ and A^ ; using the relation 
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Writing the transfer matrix as a product of sparse ma- 
trices Q, we obtained results up to N = 19 layers for 
the (11) orientation and N = 37 for the (10) orientation. 
For more details see Ref. Gal. 
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FIG. 2. Transfer-matrix results for strips of N layers with 
(10) and (11) surfaces. For the (10) and (11) orientations, 
t = 2 and 1, respectively. The middle two and lower two 
curves show results for even and odd N, respectively. 

Representative numerical data are shown in Fig. 2. 
The quantity irN^ 1 ^ extrapolates convincingly to the 
amplitudes 2, i, and 1 in Eq. (^), in agreement with the 
surface universality classes predicted above. 



The ordinary surface critical behavior expected in sys- 
tem (^|) for the (11) orientation and arbitrary values of H, 
Hi, and for the (10) orientation in the case H = H \ = 
(ideal stoichiometry) is confirmed by the transfer-matrix 
results. The corresponding data in Fig. 2 agree well with 

^ord.ord 2 

We predict extraordinary surface critical behavior for 
the (10) orientation with nonvanishing H, Hi, due to 
effective ordering surface fields, with ++ and H — orien- 
tations for antiferromagnets with odd and even N, re- 
spectively. This is also confirmed by the transfer-matrix 
results. The corresponding data in Fig. 2 are in excellent 
agreement with A ++ = \ for N odd and A' 1 = 1 for 
N even. The magnetization and order-parameter profiles 
in strips of 37 and 36 layers with (10) edges, calculated 
from the eigenvector with the largest eigenvalue , 
shown in Fig. 3 and Fig. 4, respectively. 
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FIG. 3. Magnetization profile of a strip of N — 37 layers 
with (10) surfaces and H x = 0, H = 4, K = K C (H). Sublat- 
tices A and B are indicated by □ and x , respectively. The 
order parameter (O) is defined as half the difference of the A 
and B sublattice magnetizations in adjacent layers. 




FIG. 4. Same as Fig. 3, but with N = 36. Since the surface 
spins on opposite edges belong to different sublattices, the 
order-parameter profile is antisymmetric. 

Finally we consider the implications of mean-field the- 
ory for the surface universality class. Mean-field theories 
are, of course, of interest as tractable first approxima- 
tions. More importantly, one can often infer the appro- 
priate continuum Hamiltonian 7i[(f>] for renormalization- 
group analysis from mean-field theory. The standard pro- 
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cedure is to begin with mean-field difference equations on 
a lattice and taking the continuum limit. The appropri- 
ate H[4>] is minimized by the resulting differential equa- 
tions and boundary conditions, which correspond to the 
Landau equations 8TL/ 5<f> = 0. 

As in Ref. (§, we have carried out mean-field studies, 
but with some differences in approach and conclusions. 
Here we give a brief summary. The details will be pub- 



lished separately 16 17]]. 

It turns out to be very useful to interpret the rather 
complex mean-field difference equations as a nonlinear re- 
cursive map, an approach pioneered by Pandit and Wor- 
tis For the various surface orientations considered 

above, the mapping corresponds to discrete Hamiltonian 
dynamics. The general properties of Hamiltonian flows 
have been studied extensively [ fl9|| . The nonvanishing 
order parameter profile for T > T c in the case of a sym- 
metry breaking surface orientation can be understood in 
this context. 

Our mean-field equations also provide a convenient 
starting point for the continuum approximation and, in 
our opinion, avoid some problematic aspects of earlier 
work [§. We are led, apart from irrelevant terms, to 
the familiar </> 4 Hamiltonian ([!]) for semi-infinite systems, 
with a nonzero surface field h± in the case of symmetry- 
breaking surface orientations and zero surface field oth- 
erwise. The question of surface fields in the continuum 
Hamiltonian is conveniently analyzed in the framework of 
Landau-like symmetry arguments ]2C| ] and the "method 
of concentration waves" pi] ]. As explained in Ref. Jl7| , 
the bulk amplitudes of all concentration waves that do 
not share the symmetry of the order parameter vanish 
in the high-temperature phase. In the presence of a pla- 
nar boundary translational invariance is lost, and there 
are some nonzero amplitudes. The effective surface field 
for symmetry-breaking surface orientations arises from a 
coupling of the order parameter to nonvanishing concen- 
tration waves. 

In summary, we have shown that the surface universal- 
ity class of semi-infinite binary alloys and antiferromag- 
nets depends on the surface orientation, with ordinary 
and extraordinary behavior for symmetry-preserving and 
symmetry-breaking orientations, respectively. The un- 
derlying mechanism is rather general, and Ising spins are 
not essential. We expect the orientational dependence 
to be a common feature of semi-infinite systems with (i) 
nonordering fields and (ii) second-order phase transitions 
in which the symmetry of two or more spatially distinct 
sublattices is spontaneously broken below T c in the bulk. 

We thank Henk W. J. Blote for useful correspondence. 
T. W. B. appreciates the support of the Alexander von 
Humboldt-Stiftung. This work was also supported by 
the Deutsche Forschungsgemeinschaft through Sonder- 
forschungsbereich 237. 



Present and permanent address 

M. E. Fisher, Rev. Mod. Phys. 46, 597 (1974). 

H. W. Diehl, in Phase transitions and Critical Phenomena, 

edited by C. Domb and J. L. Lebowitz (Academic, London, 

1986), Vol. 10, p. 75. 

K. Binder, in Phase Transitions and Critical Phenomena 
(Ref. §), Vol. 8, p. 1. 

J. L. Cardy, in Phase transitions and Critical Phenomena 
(Ref. §), Vol. 11, p. 55. 

Fluids with surfaces generally correspond to Ising-like 
models with surface fields. Hence the name "normal", in- 
troduced by H. W. Diehl, Phys. Rev. B 49, 2846 (1994), 
following a suggestion of M. E. Fisher. 
The order parameter and the energy density have surface 
scaling dimension x° x = i(d — 2 + rf^ x ) = d at the extraor- 
dinary transition, implying the same surface singularity 
|T — T c \ 2 ~ a as the bulk free energy and pair correlations 



that fall off as 



[7] 



[9] 



[10] 

[11] 

[12] 
[13] 

[14] 

[15] 
[16] 
[17] 
[18] 
[19] 



[20] 
[21] 



parallel to the surface. At the ordi- 



nary transition the surface scaling dimension of the energy 
density is also d, and the surface scaling dimension of the 
order parameter is an independent quantity for which no 
simple general formula is known. For more details and a 
discussion of the extraordinary-normal equivalence see T. 
W. Burkhardt and J. L. Cardy, J. Phys. A 20, 1233 (1987); 
H. W. Diehl and M. Smock, Phys. Rev. B 47, 5841 (1993); 
48, 6740(E) (1993); T. W. Burkhardt and H. W. Diehl, 
Phys. Rev. B 50, 3894 (1994); and references therein. 
For a discussion of the ordinary-normal crossover for weak 
hi see U. Ritschel and P. Czerner, Phys. Rev. Lett. 77, 
3645 (1996). 

F. Schmid, Z. Phys. B 91, 77 (1993). 

H. Dosch, Critical Phenomena at Surfaces and Interfaces, 

Springer Tracts in Modern Physics, edited by G. Holder 

and A. Niekisch (Springer, Berlin, 1992), Vol. 126. 

H. W. J. Blote and X. N. Wu, J. Phys. A 23 L627 (1990); 

X. N. Wu and F. Y. Wu, Phys. Lett. A 144, 123 (1990). 

B. M. McCoy and T. T. Wu, The Two Dimensional Ising 

Model (Harvard University Press, Cambridge, MA, 1973). 

J. L. Cardy, Nucl. Phys. B 275, 200 (1986). 

T. W. Burkhardt and T. Xue, Phys. Rev. Lett. 66, 895 

(1991). 

H. W. J. Blote, F. Y. Wu and X. N. Wu, Int. J. Mod. Phys. 
B4, 619 (1990). 

Anja Drewitz, Diplomarbeit UGH Essen (1996). 

R. Leidl, A. Drewitz, and H. W. Diehl, to be published. 

R. Leidl, to be published. 

R. Pandit and M. Wortis, Phys. Rev. B 25, 3226 (1982). 
A. J. Lichtenberg and M. A. Lieberman, Regular and 
Chaotic Dynamics, 2nd ed., Series in Applied Mathemati- 
cal Sciences, Vol. 38 (Springer, Berlin, 1992). 
L. D. Landau and E. M. Lifshitz, Statistische Physik, 
Lehrbuch der Theoretischen Physik, Vol. V, 6th ed., 
(Akademie-Verlag, Berlin, 1984), Ch. XIV. 
A. G. Khachaturyan, Prog. Mat. Sc. 22, 1 (1978). 



4 



